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Abstract 

We compute the SUSY effective hamiltonian that describes the | AS\ = 1 semilep- 
tonic decays of tau leptons. We provide analytical expressions for supersymmetric 
contribution to r — > usv T transition in mass insertion approximation. We show that 
SUSY contributions may enhance the CP asymmetry of r — > Kttv t decays by several 
orders of magnitude than the standard model expectations. However, the resulting 
asymmetry is still well below the current experimental limits obtained by CLEO 
collaborations. We emphasize that measuring CP rate asymmetry in this decay 
larger than 10 -6 would be a clear evidence of physics beyond the supersymmetric 
extensions of the standard model. 

1 Introduction 

Strangeness-changing tau lepton decays play an important role in testing the dy- 
namics of \AS\ = 1 weak interactions [1] . Determination of basic parameters of the 
Standard Model (SM) and tests of fundamental symmetries can be done using such 
tau decays. For instance, measurements of the spectral functions of tau decays into 
strange mesons have been used recently to obtain information on the mass of the 
strange quark and on the V us CKM matrix element [2]. Furthermore, searches for 
CP violation effects in the double kinematical distributions of — > Ks^v T decays 
have been performed recently by the CLEO Collaborations [3]. These exclusive de- 
cays can be used to provide further tests on the violation of the CP symmetry [3-6] . 
A 'known' CP rate asymmetry of O(10~ 3 ) has been pointed out to exist between 
t~ — > Kl s^ ^t and their CP conjugate decays [5]. On the other hand, within the 
SM, the CP rate asymmetry turns out to be negligibly small (of order 10 -12 ) in 
— ► K^ir ^ decays [6], opening a large window to consider the effects of New 
Physics contributions. 
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Supersymmetry (SUSY) is one of the most interesting candidates for physics 
beyond the SM. In SUSY models there are new sources of CP and flavor violation 
that may lead to significant impacts on the CP asymmetries of r-decays. In this 
paper we analyze SUSY contributions to the CP violating effects in r — > Kirv T 
decays. We consider the effects due to the chargino and neutralino exchanges by 
using the mass insertion approximation (MIA) which is a very effective tool for 
studying SUSY contributions to flavor changing neutral current processes (FCNC) 
in a model independent way. We take into account all the relevant operators involved 
in the effective Hamiltonian for | A5| = 1 r decays and provide analytical expression 
for the corresponding Wilson coefficients. 

The elementary process underlying \AS\ = 1 decays is r~ — ► suv T . The lowest 
order contribution to this decay in the SM is mediated by the exchange of a single 
W~ boson, which is Cabibbo suppressed. We consider in this paper the higher order 
effects induced by supersymmetry in the effective Hamiltonian to describe this low- 
energy process. The observable effects induced in some of the dominant |A5| = 1 
exclusive processes are considered. The main focus of our paper is to provide an 
specific mechanism to generate CP-violating couplings in the scalar form factor of 
r — > Kttv t decays as studied by the CLEO collaboration in [3]. 

This paper is organized as follows. In section 2 we briefly review the SM con- 
tribution to the |A5| = I r decays. As pointed out before, the resulting CP rate 
asymmetry in — > K ± tt°u t is negligible. Moreover, we show that this result re- 
mains intact in the case of minimal extension of the SM with right handed neutrinos. 
In section 3 we derive the SUSY effective hamiltonian for the t~ — > suv T transitions. 
In Section 4 we consider the effects of SUSY contributions on the branching ratios of 
two dominant exclusive |A5| = 1 decays, namely r~ — ► K~v T and r _ — ► (Kir)~v T . 
Section 5 is devoted for analyzing the SUSY contributions to the CP asymmetry 
in r — > Kirv T decay. We show that within SUSY models, one can generate the 
CP asymmetry, however it is below the experimental limits. Finally, we give our 
main conclusions in section 6. We have also included an Appendix to provide the 
complete expressions of the Wilson coefficients derived from SUSY. 



2 | AS | = 1 r decays in the Standard Model 

Strangeness-changing \AS\ = 1 decays of tau leptons are driven by the r~ — > usv T 
elementary process. In the SM they occur at the tree-level, as shown in Fig. 1. The 
SM effective Hamiltonian underlying these decays is given by 

G F 

T-Lsm = -j=V us (v T 'y fl LT)(s'y fl Lu) , (I) 

where V us is the its Cabibbo-Kobayashi-Maskawa (CKM) matrix element and L,R = 
1 T 75- The amplitudes for the dominant exclusive processes derived from this 



2 



Figure 1: SM tree- level contributions to r — > usv T transition. 



Hamiltonian are: 

A SM (T-(p)^K-(q)u T (p')) 

A SM (T-(p)^K{q)irtf)v T tf) 



■i—^V us f K m T r 

i^v r 

1 Vus^K 



+ rn T f s (t) r 



f v {t) Hp)^Lt{p) 



(3) 



where letters within parenthesis denote the momenta of the particles, fx is the K~ 
decay constant, t = {q + q') 2 is the square of the momentum transfer, C K = l(l/v2) 
for K°7r~(K~7r°) state, r = v(p')Rt(p), A 2 = m\ - m 2 , and 



Qtx = (q 



A 2 



(? + </: 



(4) 



In the SM, the two-body decay of Eq.(2) is a clean prediction if one uses fx — 
159.8 MeV from K~ — > pTVn decay [8]. Since New Physics can affect r — > Kv 
and K — > [iv decays in a non-universal way, these decays can be used to obtain 
interesting bounds on new physics couplings. 

On the other hand, the three-body decay of Eq.(3) can exhibit eventually the 
effects of CP violation [3-6]. However, the decay rate of this process is given by [7] 



r (r -> KlTVr) 



Glm 



T l Vus\ 2 Ism, 



768tt 3 



(5) 



where 



Ism 



1 



dt 



+ SI/^A 4 ^,™ 2 ,,™ 2 )] 1 / 2 



\fv\ 2 { 



1 + 



2/ 



mi 



2^3/2 



(6) 



The function X(x, y, z) is given by X(x, y, z) = x 2 +y 2 +z 2 —2xy—2xz—2yz. It is clear 
from the expression of T(t — ► Kirv T ) that within the SM, the direct CP asymmetry 
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identically vanishes. As is well known, a necessary condition to generate a CP 
rate asymmetry is that at least two terms of the amplitude for a given physical 
process have different weak and strong phases. However, in the SM the relative 
weak phase between the scalar fg and the vector fy form factors of r — ► Kttu is 
zero. Furthermore, since the form factors fv,s(t) m Eq. (3) belong respectively 
to the (orthogonal) I = 1 and / = angular momentum configurations of the Kir 
system, the fgfv term in the squared amplitude vanish upon the integration over 
the variable u = (p — p r ) 2 , therefore the CP- violating terms vanish in the integrated 
rate T and in the hadronic spectrum dT/dt. Thus, within the SM CP violating 
effects in this three-body channel can manifest only in the double differential decay 
distribution where the interference of fy and fg is present (see ref. [3] ) . 

A different mechanism to generate a CP rate asymmetry in the SM for — ► 
K^ir^Vr was considered in [6]. In this case the two amplitudes with different weak 
and strong phases contribute to the same I = 1 angular momentum configuration. 
This asymmetry turns out to be negligibly small since it is suppressed by the CKM 
factor Vu — 1CT 3 and also by a higher order suppresion factor g 2 /4irMyy ~ 10~ 8 . 
Thus, the resulting CP rate asymmetry is expected to be negligible, as confirmed in 
Ref. [6]. Therefore, this decay can be suitable to search for the effects of CP violation 
induced by New Physics. 

The minimal extension of the SM with right handed neutrinos z^SM, where non 
vanishing neutrino masses can be obtained, allows for a new source of the CP 
violation through the Umns mixing matrix. In this scenario, the amplitude of the 
decay r~ — > K~ir°v T will be given by 

A vSM {t~ -+ K~tt°v t ) = (U* MNS ) 33 \A SM \ e l5sM . (7) 

It is remarkable that, although the amplitude J[ uS - M has a weak CP violating phase, 
the CP asymmetry still vanishes. Therefore, any measurement of a non-vanishing 
CP asymmetry will be a hint for a new physics beyond the SM. In the rest of the 
paper we will focus on the NP contributions to CP violation in the three-body decay 
induced by SUSY. 



3 SUSY contributions to \AS\ = 1 r decays 

The effective Hamiltonian H e ff derived from SUSY can be expressed as 

H eff = ( ^=V us '£Ci(riQi(Li), (8) 

where Cj are the Wilson coefficients and Qi are the relevant local operators at low 
energy scale fx ~ m T . The operators are given by 

Qi = {v^Lt){s 1ix Lu), (9) 
Q 2 = (v^Lt)(s^Ru), (10) 
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Q 3 = (uRt)(sLu), (11) 
g 4 = (PRt){sRu), (12) 
Q 5 = ipa^Rr^sa^Ru). (13) 

where L,R are as defined in the previous section and = |[7 /i ,7 I/ ]. The Wilson 
coefficients Q, at the electroweak scale, can be expressed as C{ = Cf M + Cf u 
where Cf M are given by 

Cf M = 1, 

C 2 5 if 4l5 = 0. (14) 

SUSY contributions to the Hamiltonian of t~ — ► usv T transitions can be generated 
through two topological box diagrams as shown in Figs. (2,3). Other SUSY contri- 
butions (vertex corrections) are suppressed either due to small Yukawa couplings of 
light quarks or because they have the same structure as the SM in the hadronic ver- 
tex. In our computations of Wilson coefficients, we will work in the mass insertion 




Figure 2: SUSY box contributions to r — ► usv T transition. 

approximation (MIA), where gluino and neutralino are flavor diagonal. Denoting by 
{/S. AB ) a i, the Ojff-diagonal terms in the sfermion mass matrices where A,B indicate 
chirality, A,B = (L, R), the A — B sfermion propagator can be expanded as 

(f%%) = itfl - mH - A/,))- * + + Q(A% (15) 

where / denotes any scalar fermion, a,b = (1,2,3) are flavor indices, I is the unit 
matrix, and rh is the average sfermion mass. It is convenient to define a dimension- 
less quantity (5 AB ) a b = (A^ B ) a h/m 2 . As long as (A^ B ) a b is smaller than fh? we can 
consider only the first order term in (5 AB ) a b of the sfermion propagator expansion. 
In our analysis we will keep only terms proportional to the third generation Yukawa 
couplings and terms of order A where A = V us . 
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Figure 3: Crossed diagrams of Figure 2. 



The complete expressions for the Wilson coefficients Cj at mw scale induced by 
SUSY computed from Figs. (2,3) can be found in the Appendix. As can be seen from 
this appendix, the Cj are given in terms of several mass insertions that represent 
the flavor transitions between different generations of quarks or leptons. In general, 
these mass insertions are complex and of order one. However, the experimental 
limits of several flavor changing neutral currents impose severe constraints on most of 
these mass insertions. In the following, we summarize all the important constraints 
on the relevant mass insertions for our process. 

1. From the experimental measurements of BR(/j, — > cy) < 1.2 x 10~ , the 
following bounds on |(<^ 2 )ab| and | (<5^ 2 )as I are obtained [10]: For M\ ~ M 2 = 



2. From BR(t — > ^7) < 1.1 x 10 6 , one gets the following constraint on |(5 2 3)i,_r)| 



100 GeV and fx = m { = 200 GeV, 



\(S[ 2 )ll\ < 10- 3 , 
1(^)^1 <6x 10 



\(S{ 2 )lr\<10- 6 , (16) 
M 3 ) LL \ < 4 x 10- 4 , \(S^) LL \ < 7 x 10- 4 . (17) 



[13]: 



and from BR{t 



l(4M £2 x 10~ 2 , 
ei) < 2.7 x 10" 6 , one finds [13]: 



(18) 



IK 3 M<ixio 



1 



(19) 






(21) 
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4. The mass insertion (<5" 2 )ab ar e constrained by the AMo as follows [14]: 

\(5i2)ll\ < 1.7 x 1(T 2 , m 2 ) LR \ < 2.4 x 1(T 2 . (22) 

Here, three comments are in order, i) Due to the Hermiticity of the LL sector 
in the sfermion mass matrix, (5 ab )ll = ($ab)ll = (^ba)lL' wnere A,B = 1,2,3. 
ii) The above constraints imposed on the mass insertions (5 ab )ll,lr are derived 
from supersymmetric contributions through exchange of gluino or neutralino which 
preserves chirality, therefore same constraints are also imposed on the mass inser- 
tions {5 a ^ b )rr^rl- Hi) The mass insertions [8 ab )lr(rl) are n °t, in general, related 
to the mass insertions (S ba )lr(rl)- Taking the above constraints into account, one 
finds that the dominant contribution to the r~ — > usu T is given in terms of (5% 2 )lr, 
(^32) RLi {&2i)RL anci {&2i)lr- Notice that the effective Hamiltonian (eq. 8) derived 
in this section can induce supersymmetric effects in all the |Ai>| = 1 exclusive r 
lepton decay. In the following section we consider two examples. 

4 Constraints from two-body r-decays 

In this section we analyze the possible constraints that may be imposed on the 
SUSY contributions to r — > s u v T from the exclusive |A5| = 1 decay: T~(p) — > 
K~(q) v T (p'). The decay amplitude considering effects of SUSY contributions reads: 

A(t~ -»• K~v T ) = Asm + Asusy , (23) 

which can be written explicitly as: 

A ( T - _> k-u t ) = i^fusfK m T (is(p')RT(p)) {1 + 5 S usy (r)} . (24) 

The decay K~ — > /it _ z/, which fixes in the absence of new physics, would also be 
modified by the effects of new physics: 

A(K~ - ^-i/) = i^|u us / x m M (u(p')Rfi(p)) {1 + 5 5[ /5y(^)} • (25) 

In order to estimate the size of SUSY contributions in such decays, one defines the 
ratio: 

r(r -> KvtW) 

r(if-H7)) 

2m K m^ / m 2\ 



A' , 



In the above equation (7) means that complete SM radiative corrections of 0(a) 
have been included. The long-distance radiative corrections, which do not cancel 
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in the above ratio, were computed in [9] and reads SR T / K = (0.90l ' 26 )% . Using 
the experimental rates for the involved decays: T(r — > Kv) = (2.36 ± 0.08) x 
10 10 s~ 1 ,r(K -► fiv) = (0.5118 ±0.0018) x lO^ -1 [8] into eq. (21) we get: 

Re[5 SU SY{K) - S S usy(t)] = 0.02 ± 0.03. (26) 

This equation is actually a model- independent result. Using the expression for the 
Wilson coefficient induced by SUSY, it is easy to see that the dominant contribution 
to Re{5susY{T~) — Ssusy(K)) comes from Cf u . In terms of the SUSY parameters, 
it can be translated as follows 

Re[S SU SY(K) - 5 SU sy(t)] ^ Re[Cf USY {K) - Cf USY {r)\ (27) 

where Cf USY (K,t) are respectively the Wilson coefficients corresponding to the 
operators responsible of K — > fii/n and r — > Knv T . Using as input parameters 
Mi = 100, M 2 = 200 GeV and p = M- q = 400 GeV and tan (3 ~ 20, one gets 

Re[Cf USY {K) - Cf USY {r)} ~ -0.03(<5 2 U 1 ) LL - 0.006(4) LL (28) 

where we keep the dominant contributions only. Using the bounds on the <5's given 
in previous section, the terms which is mostly unconstraint is (S 1 23 )ll- It is clear 
that still we are far from the experimental limit given in (26) but it is important to 
notice that experimental data should strongly improve in a close future[12] . 



5 SUSY contribution to CP asymmetry in r — > Ktiv t 

Having analyzed the constraints from two-body r-decays imposed on the supersym- 
metric contributions to r — > s u v T transition, now we can study the supersymmetric 
effects on the CP violation in the three-body decay r — ► Kttu t . We will show that 
although supersymmetry enhances the asymmetry of this process by many order of 
magnitude than the SM expectation, the resulting CP asymmetries are still smaller 
than the current experimental reaches. 

Given the spin-parity properties of the Kir system, we can write the total am- 
plitude (SM and SUSY) of the r(p) -► K(q)ir(q')u T (p') decay as 

•At(t — ► Kttu) = ^^f(l + Ci)(ir7r|s 7 ^|0)P(p / )7' i ^(p) (29) 

V2 L 

+ (C 3 + Ci)(K7r\su\0)i?(p')RT(p) + C 5 (K7r\sa^u\0)i>(p')a^ Rr(p) 

where Cj stand for Cf USY , since the Cf M are explicitly included . It is now clear 
that the resulting CP asymmetry depends on the relative ratio among the SUSY 
Wilson coefficients. For example, in case that C\ is giving the dominant SUSY con- 
tribution and C3 4 5 effects can be negligible, then the CP asymmetry of r — > Kixv T 
will vanish identically as in the SM. We consider the following two interesting sce- 
narios: 
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i) The case of C3 or C4 gives relevant contributions while C5 is negligible. In 
this case, SUSY induces a relative weak phase between the vector and scalar form 
factors describing this process. 

it) The case of C5 gives relevant contributions while 63,4 are negligible. In this 
case, SUSY induces a relative weak phase between the vector and tensor form fac- 
tors. 



Let us start by analyzing the CP asymmetry in the first scenario. Using the 
definition of the hadronic matrix element introduced in Eq.(3): 

(Kir\s llx u\0) = C K {f v {t)Q^ + f s (t)(q + , (30) 

we can obtain the hadronic matrix element of the scalar current by taking the 
divergence in the usual form: 

(Kn\su\0) = ° Kt f s (t) , (31) 
m s - m u 

where vn, s u denote s, u current quark masses. Thus, we finally get the amplitude: 



v2 

f v Q»u(jp')^L U (jp) + 



m T + 



C 3 + C 4 \ t 



1 + C\ ) m s — m t 



f s u(p')Ru(p) 



It is remarkable that in this case, SUSY effects just modify the normalization of the 
SM amplitude and the relative size of the vector and scalar contributions. 

When we compare this expression with the decay amplitude given in Eq. (2) of 
Ref. [3]: 

A(r~ - Kirvr) ~ u{p')l»Lu(p)f v Q» + Au( P ')Ru(p)f s M , (32) 
where M = 1 GeV is a normalization mass scale, we obtain the relation 

The first term in the r.h.s of this equation is the usual contribution of the SM, 
which is real, and the second term arises from the SUSY contributions and contains 
a CP-violating phase, hence it can generate a non-vanishing CP asymmetry. 
The square of the matrix element becomes: 

^|^l| 2 ~ \f v \ 2 (2 P .Qp'.Q-p.p'Q 2 ) + \A\ 2 \fs(t)\ 2 M 2 p.p' 

pols 

+2ReA • Re(f s fv)Mm T p'.Q - 2/mA • Im(f s fv)Mm T p'.Q (34) 

The last term in the previous equation is odd under a CP transformation but we 
should notice that the last two terms disappear once we integrate on the kinematical 
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variable u of the process in consideration. This means that it is not possible to 
generate a CP asymmetry in total decay rates corresponding to this process. So 
the only way to generate a CP asymmetry is to look for the double differential 
distribution (d 2 T / dudt) or a variant of it as CLEO collaboration did in Ref.[3]. The 
CLEO collaboration has recently studied the ratio of CP-odd to CP-even terms of 
this squared amplitude for — ► Ks^v T decays and has obtained the following 
bound: -0.172 < Im(A) < 0.067 at 90% C.L. Using Eq.(33), we can translate this 
bound into: 

-0.010 < Im ^ i + c A ) ^ °- 004 > ( 35 ) 

where we have used m s —m u = 100 MeV, and the average value (t) ~ (1332.8 MeV) 2 . 
Now, as input parameters Mi = 100 and M 2 = 200 GeV and fj, = Mg = 400 GeV 
and tan (3 = 20, one gets 



Im I 



- 1-3 x W- 5 Im(5 d 21 ) RL (36) 



Still the experimental bound is too loose to give us information on (S^rl- Notice 
however that forthcoming measurements of the CP asymmetry in the Kstt channel 
will be significantly improved at B factories [12] since their data sample is larger by 
two order of magnitude than the one used by CLEO [3] in their analysis. 

Now, we turn to take into consideration the O5 operator which is naturally in- 
duced by SUSY corrections to Wilson coefficients. This operator could interfere with 
0\ operator which contains SM contributions and the strong phases. So in princi- 
ple, using this interference between O5 and 0\, it should be possible to generate a 
CP asymmetry directly in the total decay rate which is completely forbidden in SM. 
Let us keep the dominant contribution to r~ — > (Kn)~v T and the O5 contribution: 



+ 



j^{Kir\s*p V u\0)u(pr)<T^Ru(p)} (37) 



In this expression, we have neglected fs effect since its contribution to total decay 
rate is numerically small (around 3% at most, see [7]) and conserves CP. The most 
general form of the antisymmetric matrix element of the hadronic tensor current is 
given by 

in 

(Kir\sa^ v u\0) = [(p^imY ~ {p*Y{pkT\ (38) 

m K 

where a is a dimensionless quantity which fixes the scale of the hadronic matrix 
element. It is important to remember that fv(t) contains the strong phases as it 
can be parametrized [7]: 

_ /v(»K- 

m K* —t — irriK**- K* 
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The tensor form factor which is given by 

A-irk <40) 

has no strong phases but of course could have a weak phase (arises either from C\ 
or C5) and it can be at most a slightly varying function of t. So, one can compute 
now the CP asymmetry in total decay rate: 

r(r~ K~tt°u t ) - r(r+ -► K + ir°v T ) 
aCP ~ T(t- FA^ + T(r+ i^+vrV) ( ^ 

x {m 2 (i + A 2 ) 2 + m 2 A + [(t - m 2 ) 2 - m^(i + m 2 )] (/ - A 2 )} (42) 
where A, A and Tsm are given in section 2. Integrating numerically on t, one gets 

a C p ^ ^ Jm ^ ( 43 ) 
~ 1.4 x 10 _7 a M^Ilb (44) 

where to get the last equation, we use the same SUSY parameters as before. Again 
within SUSY extensions of the SM, this CP asymmetry is small (even if it is practi- 
cally five orders of magnitude bigger than the one expected in SM[6]) . Clearly the 
observation of a CP asymmetry in this channel at a range bigger than 1CT 6 will be 
not only a clear evidence of Physics beyond Standard Model but also an evidence 
we need Physics beyond Supersymmetric extensions of Standard Model. 



6 Conclusion 

In summary, in this paper we have computed the effective hamiltonian derived from 
SUSY for \AS\ = 1 tau lepton decays using the mass insertion approximation. Al- 
though experimental data for such decays are not precise enough at present to give 
constraints on the fundamental parameters of SUSY, we have shown how physics 
beyond standard model as supersymmetrix extensions of the SM could induce CP 
violating asymmetry in the double differential distribution as CLEO collaboration 
did in ref.[3] and could also induce CP asymmetry in total decay rate due to inter- 
ference between O5 and 0\ operators. We have argued that any CP asymmetry in 
the channel under consideration bigger than 10 -6 will be a clear evidence of not only 
Physics beyond Standard Model but also an evidence of Physics beyond SUSY ex- 
tensions of the SM. We also provided model-independent constraint on New Physics 
contribution to r — > Kv. In particular, it is interesting to observe that SUSY can 
provide a specific mechanism to generate a CP-violating term in the probability 
distribution of r — > Ki:v T decays. Forthcoming and more precise data for observ- 
ables of the exclusive processes considered in this paper will either provide better 
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constraints on SUSY parameters or give a mechanism to explain discrepancies with 
the SM if they are observed. 
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8 Appendix 

SUSY contributions to Wilson coefficients of r — > suv T 

Here we provide the complete expressions for the supersymmeric contributions, at 
leading order in MIA, for the Wilson coefficients of r~ —> suv T transition, Ci(Mw), 
i = 1, .., 5. As mentioned in section 3, the dominant SUSY contributions are given 
by chargino-neutralino box diagram exchang illustrated in Fig. 2. 
The effective Hamiltonian -He// derived from SUSY can be expressed as 

^V^di^Qiifj), (45) 



H eff 



where Cj are the dimensionless Wilson coefficients and Qi are the relevant local 
operators at low energy scale fi ~ m T . The operators are given by 

(47) 
(48) 
(49) 
(50) 
(51) 

In terms of the vertex, one can write the complete vertex as a product of the 
vertex coming from leptonic sector and of the vertex coming from hadronic sector. 
In this respect we can also write the Wilson coefficients as 

^ = C kr-X-) ( C i(r- X --s) +C !((r- X --u)) 
+ C kr-X°) ( C i{T-x°s) + C i((T- X °-u) 

12 



Qi = 




Q2 = 


(i/7 M Lr)(s7 /i J?u), 


Qs = 


(uRt)(sLu), 


Q 4 = 


{uRt){sRu), 


Q 5 = 


{vo^Rt^sg^Ru) 



where the C\ is due to the leptonic vertex and Cf is from the quark sectors. If we 
expand C/ 9 in terms of the mass insertions, one finds that the leading contributions 
are given by 

+ ^l(r-X _ ) ipl{T- X --s)^ n ^ Xi ' X j) + ^l{T- X --u) In ^ Xi -' X i)) 
+ C^_ x0 ^ (cf^_ x0 _^I n {Xi, Xj) + C^}_ x0 _ u ^I n {Xi, Xj^j 

+ ^l(-r-X-) (^'l(r-x--s)^ n ^' X J^ + ^(r-x"-")^ 1 ^*' 

+ Ci( T l_ x O) (Ci( T - x O_ s )-^n(^i) + Cl( r - x 0-u)^n( x i5 (52) 

+ 0(5 2 ). 



With 



= ^ (JV*2 + tan^JVji) ^(^5)33 - y 2 T N* 3 U* 2 (U* MNS ) 33 , (53) 
^it°- x o- u ) = (^(^2 + ^tan^iV*)^^)^, (54) 
Ci( ( T °-x°- S ) = (^)^(^ 2 -^tan^iV a )c/; i (VS i , M ) 12 , (55) 

2 

C^o) = -^=(^ + tan^iVA)^l(^M7V5)a3(^L)a3 (56) 

+^ W 2 + tan^TV* ) ^(^5)33(^)33(^)33 

+g(heh3N* 3 U* 1 (5 l LR ) a3 (U* MNS ) a3 
-(h e )l 3 N* 3 U* 2 (S l RR )3 3 (Ul INS ) 3 3, 



and 



Cit- x o_ s) = (^)(^(^ 2 -^tan^iv, 1 )c/; 1 (y^ M ) la (^ L )2a 



+JL (N i2 - Itan^iV a j ^2(^^)13(^)33(^)23] , (57) 

<?it- x o_ u) = (^)^W2 + ^taa^ 1 )(^ Af ) o2 («^ £ ) ol ^ 1 

+ lt a ne w N* 1 )(h u ) 33 V j2 (VS KM h2(5 u RL h^ . (58) 



CiCr-x-) = 7^^1(^5)33(^2 - tan^iV a ) - (K)\ 3 V* 2 (U* MNS ) 33 N i4 , (59) 
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9(0) _ / 1 n / 5 2 , 



l((r- X — ) = (g)|-^(^2-3^^l)(^M)l2l/;il, (61) 



,/(!) 5 2 T , , , 



q^_ x _j = ^=F/ 1 (^ JV5 ) 3 a(iV i 2-tan^iV a )(^)3 a 

+gV* 1 (h I/ ) 33 N i4 {6 RL ) 3a (U* MNS ) 3a 

~~^( Ni2 ~ tan0 ™ N il) V M h ") aa ( U MNsha(5LR)3a 
-(KhsN^iKUitRRhaiUhNsha. (62) 

In case of decoupling of the sneutrino-right, the last terms are strongly suppressed. 



Cf$- X - U ) = (^)(^^i(^2 + ^tan^)(^ M ) a2 (^ L ) al 

-^=(h u ) 33 V j2 (VS KM ) 32 (5 u RL ) 31 (N* 2 + itan^AT*) ] (0:5) 



U l(r-x"- S ) ~ ^> 



5 , Ar 1 



§) I (JVi2 - §tan^iV il )(^ M ) la C/; i (<5i L ) 2a 



+^(^<2 " 3 tan^iV il )(^Af)i3^ 2 (M33&)23j • (64) 
The contribution to C 2 is found to vanish identically, i.e., 

C 2 = 0, (65) 

A, _ WW r ,V T .\ , WW T (t- t ) 

° 3 ~~ 3(r-x") 3(t-x"-s) inl > X *' L/ 3(r-x ) O 3(T-x°-s)-''' t V X «' X J^ 

+0(5 2 ), (66) 
where I n (xi,Xj) is defined below and Xi = rn^ ± /rh 2 and Xj = m^ /rh 2 . 

C$_ x0) = 5(^)33^3^*1(^5)33 (67) 
-gV2tane w N il U* 2 (h e ) 33 (U* MNS ) 33 , (68) 

Cjjl^) = -(h e ) 33 ^=(N i2 - t a n9 w N a )U ]2 (U* MNS ) 33 , (69) 



°3(r- X -- S ) " i oJ 



^^tan^iV*^ 2 (/ l(i )33(^ M )i3(4i ? )23) , (70) 
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-^tan^£/* 2 AT* (V^M^rMMss) • (71) 



C4 — C 4 ( T _ x _- ) C| ( - T _ x __ u - ) 7 n (xj, Xj) + ^4( T _ x O)C'4(( T _ x o_ u )-fn(^i) »j) + 0{5 2 ), (72) 
where I n (xi,Xj) is given below. 



C? ( ° T ) _ x - ) = -(^)33^(^2-tan^iV a )^ 2 (C/^ Ar5 )33 

= C$- x -y (73) 



fd^) = 5(^)33^3^(^5)33 

-gV2tane w N il U* 2 (h e ) 33 (U* MNS ) 33 , (74) 



^- x o_ u) = (-l)^-l^tan9 w N il V jl (VS KM )a2(S u LR ) a i 

+^ ta ne w N il V j2 (h u ) 33 (VS KM )32(5 u RR hi) , (75) 



= (-^)(-^ tan ^^l^l(^M)a2&)al 

+ i^ tan ^y. 2 7V a (^)33(^ M )3 2 (^) 31 ) . (76) 

^5(t-x°-«) = "7^4(t-x°-«)' (^7) 



1 * 

4 

The loop integrals I n (xi,Xj) and I n (xi,Xj) are defined as follows: 



C '5(t-x--m) - 7 C 4H"-«)' (78) 



w x _ 1 / 1 \ / A -Xj- xflogxj , 



~ ^xjx] ( 1 \ ( x\ -Xj- Xjlogxj \ 

Hxi, Xj ) - lg7r2 ^ 2 y x ._ x .j y (1-x,) 2 W" x j))> 
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(Xi <-> Xj)) , 



In(%ii Xj) 




+ 3xi + 2xilogxi 




(79) 
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